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Abstract. In this paper, we establish a novel approach to proving existence of non-negative 
weak solutions for degenerate parabolic equations of fourth order, like the Cahn-Hilliard and 
certain thin film equations. The considered evolution equations are in the form of a gradient 
flow for a perturbed Dirichlct energy with respect to a Wasserstein-like transport metric, and 
weak solutions are obtained as curves of maximal slope. Our main assumption is that the 
mobility of the particles is a concave function of their spatial density. A qualitative difference 
of our approach to previous ones is that essential properties of the solution — non-negativity, 
conservation of the total mass and dissipation of the energy — are automatically guaranteed by 
the construction from minimizing movements in the energy landscape. 



1. Introduction and statement of main results 

This paper is concerned with the following class of initial-boundary value problems for non- 
linear fourth order parabolic equations, 

d t u = -div(m(u)D(Au-G'(u))) in(0,oo)xft, (1) 

n-Du = 0, n- (m(u)D(Au-G'(u))) =0 on(0,oo)x9O, (2) 

u(0,x) = uq(x), in fl. (3) 

The problem (UJ)-© is posed on a bounded, smooth convex domain 51 C R d . n denotes the 
normal vector field to the boundary d£l. The sought solution u : [0, oo) x 17 — > R is subject to the 
constraint < u(t, x) < M, where either M > is a given number, or M = +oo. The mobility m 
is a non- negative concave function m : (0, M) — > R + that vanishes at 0, and also at M if M < oo. 
The mobility m and the free energy G : (0, M) — > R are subject to certain regularity assumptions, 
specified below. Introducing the pressure P satisfying 

P'{s) = m( S )G"( S ), (4) 

equation ((TJ) can be rewritten in the more familiar form 

d t u = - div (m(u) D Aw) + AP(u). (5) 

Equations of the form ([T )) or ([5 jl arise, for instance, as hydrodynamic approximation to models for 
many-particle systems in gas dynamics, and also in lubrication theory. In particular, the classical 
Cahn-Hilliard equation for phase separation in a binary alloy as well as the (de)stabilized thin 
film equation are of the shape (JTJ); we comment on these special cases further below. The value 
of the solution u(t, x) represents a particle density, or the fraction of one component of a binary 
alloy (in the case of the Cahn-Hilliard equation), or the height of the film (in the case of the thin 
film equation) at time t > and location x G fl. 

There is a rich literature on the mathematical structure of Cahn-Hilliard, thin film and related 
equations. In particular, the techniques developed in the seminal papers by Elliott and Garcke 
[14] and by Bernis and Friedman [2] have been proven extremely powerful to carry out existence 
analysis, and have been extended by many other authors afterwards. As a core feature, these 
techniques allow to replace ((TJ by a family of regularized problems with smooth solutions ug that 
satisfy certain bounds which produce the desired constraint < u < M in the limit S J. 0. We 
emphasize that this bound does not come for free since solutions to fourth order equations do not 
obey comparison principles in general. 
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More specifically, in the existence proof for the Cahn-Hilliard equation [14], the degenerate 
mobility m is replaced by a strictly positive approximation mj defined on all R. The resulting 
parabolic problems are non-degenerate and possess global and smooth solutions us, which, how- 
ever, may attain arbitrary real values. Using additional a priori estimates, it is then shown that 
the integral of us in the region where u < or u > M converges to zero as mj approaches m, 
which yields < u(t) < M in the limit. In fact, a corollary of this method of proof is that 
solutions to (JIJ with a sufficiently degenerate mobility function m preserve the strict inequalities 
< u(t) < M for all times t > 0. This property has been used in the existence proofs for thin 
film equations [5] [3] , where the original mobility is approximated by very degenerate mj . 

The techniques from [2J [T3] rely on the dissipation of certain Lyapunov functionals by solutions 
to fl}. One distinguished Lyapunov functional is a perturbed Dirichlet energy, which is defined 
on functions u G ii^ 1 (SI) with < u < M by 

e M = t; / |Dw(x)| 2 dx+ f G(u(x))dx. (6) 
2 Jo Jn 

This energy and its dissipation provide regularity estimates. Another Lyapunov functional intro- 
duced in [2J [H] and since then widely used the literature is 

Ukt] := { U(u(x))dx with U"(s) = — (7) 
Jn m (s) 

If the mobility m(s) degenerates sufficiently strongly for s i and s f M, then U allows to control 
the solution u in the zones where u is close to or M . The functional U has thus become a key 
tool for proving the bounds < u < M. 

Here we develop an alternative approach to existence which avoids the cumbersome discussion 
of the propagation of the bound < u < M and shows a new interesting variational structure 
behind equations of the form ([lj . Our starting point is the classical observation that (TT|) is in the 
shape of a gradient flow for E on the space of non-negative density functions of fixed mass. On a 
purely formal level, the corresponding metric tensor is readily determined: to a tangential vector 
v := d s p(0) to a smooth curve p : (— e,e) — > L 1 (S1) of strictly positive densities p(s) at po — p(0), 
it assigns the length 

||v|| 2 = / | D ip(x)\ 2 m(po(x))dx, with <p satisfying — div (m(po(x)) D ip(x))) = v in fl (8) 
Jn 

and variational boundary conditions on dCl. 

In the particular case of a constant mobility m = 1, this is simply the dual of the Sobolev 
seminorm in W ' (fl), and one can work in the the well know setting of gradient flows in Hilbert 
spaces, see e.g. [3J. For the linear mobility m(s) = s, the tensor is induced by a non-Hilbertian 
metric, namely the celebrated L 2 -Wasserstein distance, see [26]. In this framework, weak solutions 
to specific cases of JT]) have been obtained as curves of steepest descent in the energy landscape 
of E; see [T71H1] for respective results on the Hele-Shaw flow. 

For more general mobilities, the existence proof presented here seems to be the first based on 
the gradient flow structure of ([T]) with respect to a metric that is not the L 2 -Wasserstein distance 
nor a flat Hilbertian one. It has been proven only recently by Dolbeault, Nazaret and the third 
author [12] that even for certain non-linear mobilities m, the formal metric structure indicated in 
([5]) still leads to a genuine metric W m on the space of positive measures. One needs to assume, 
however, that m is a concave function to get nice analytic and geometric properties of W m : they 
have been studied in [7] and [53] and we review selected results in Section [2] below. 

The goal of this paper is to prove rigorously that weak solutions to ([I])-© can be obtained 
by the variational minimizing movement/ JKO scheme under suitable conditions on the nonlinear 
concave function m. The terminology minimizing movement scheme is due to De Giorgi [10,, 
whereas JKO scheme enters in common use after the paper [19j . Preservation of the total mass, 
dissipation of the energy and, most notably, non- negativity of the density along the solution are 
direct consequences of the applied construction: the solution is a weak limit of time-discrete 
energy minimizing curves that lie in the convex cone of non- negative densities. The difficulty of 
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this approach consists in proving a posteriori that the curve of maximal slope indeed corresponds 
to a weak solution. For this, a priori estimates resulting from the dissipation of U are employed. 

1.1. Hypotheses. We recall that 

either M > is a given number, orM = +oo. 

All appearing measures are assumed to be absolutely continuous with respect to the Lebesgue 
measure J? d , and we identify them with their Lebesgue densities on fi. The densities have fixed 
total mass m > 0, and are bounded from above by M if the latter is finite. Thus our ambient 
space will be the metric space (Jf(f2),W m ) where 

X(Q) := ju e L 1 ^) : < u < M a.e. in Q, ^nda; = m}. (9) 

The possible choices for mobility functions m : (0, M) —> M.+ := (0, oo) are subject to the following 
conditions: 



m is concave, m e G°°(0, M), m > in (0, M), 

m(0) := limm(s) = 0, m(Af) := lim m(s) = if M < +oo. 

sj.0 sfM 



(M) 



Moreover, we say that the mobility m satisfies a LSC condition (i.e. m is Lipschitz and m 2 is 
Semi-Convex) if 

sup |m'(s)| < +00 and sup ( - m"(s)m(s)) < +oo. (M-LSC) 

se(0,M) sS(0,M) 

The restriction to concave mobilities in (|M| is necessary, since only for those, the corresponding 
metric W m is well-defined. Notice that this hypothesis is somewhat opposite to the one made 
in [14] . where an asymptotic behaviour m(s) ~ s a for s —¥ with a > 1 has been assumed. 
Typical examples for mobility functions with finite M > are m(r) = r(M — r), or, more 
generally, m(r) = r a °(M — r) ai with exponents ao, a\ € (0, 1]. These mobilities satisfy (|M-LSC|) 
iff ao = ot\ = 1. In the case that M = +oo, the mobility m is nondecreasing because it is concave 
and strictly positive in (0, +00). Typical examples are m(r) = r a with a € (0, 1]; such a mobility 
is LSC only in the Wasserstein case a = 1. 

Concerning the free energy G 6 C 2 (0, M) and the associated pressure P with we assume 
that there exist a constant C > and an exponent q > 2 with q < 2d(d — 4) if d > 4 such that 

mG" > -C in (0,M), PeC°([0,M]) if M < oo, 

PCsl (G) 
mG" > -C(l + m) in (0,oo), P G G°([0, 00)), lim ^— - = if M = +00. 

s->oo si + |G(s)| 

The condition ([G} yields in particular (see ij2.4l) 

G e G°([0, M]) if M < 00; G(s) > -G(l + s 2 ) for every s > if M = +00. (10) 

Examples for sensible choices of G (and P) fulfilling these assumptions are given after the state- 
ments of our main results. 

1.2. The minimizing movement approximation and the existence result for LSC mobil- 
ities. The minimizing movement/ JKO scheme is a variational algorithm to obtain a time-discrete 
approximation (of given step size r > 0) to a curve of steepest descent, see pQ- In the situation at 
hand, we start from the initial condition uq € X(il) with E[ito] < +00 and define inductively 

u° T := uo, < +1 := argmin*; 1 G X(fi) where *; l (v) := ^-W m «,w) 2 + E[«], (11) 

2r 

and we set E[u] := +00 if u ^ H 1 ^). The approximation u r : [0, 00) — >• X(fi) is defined by 
constant interpolation, using u T (t) — u™ for [n — 1)t < t < nr. 

Theorem 1. Assume that Vl is a smooth bounded convex open subset ofM. d , the mobility function 
m satisfies (fM|) and (jM-LSCp . and the free energy G satisfies ([Gj) . 
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Then, for any initial condition uq G X(fi) of finite energy E[uo] < +00, the scheme (jTTJ) admits 
time-discrete solutions u T for all r > 0. For every sequence r„ 4 there exists a subsequence, still 
denoted by r n , and a function u satisfying: 

u G Lf oc ([0, 00); H 2 (Q)) n C° ([0, 00); H\Q)) n ACL([0, 00); X(it)) (12) 

u Tft (t) — > u(t) strongly in L 2 (f2) and weakly in H (Q) for all t G [0, +00), (13) 

u Tn -> u strong in L 2 (0, T; H 1 ^)) and weakly in L 2 {0, T; H 2 {n)) for all T > 0. (14) 

The energy satisfies the bound 

E[u(t)] < E[u ] for all t > 0, (15) 
there exists a decreasing function (p : [0, +00) — > R such that 

<p(t) > E[u(t)] for all t > 0, (16) 

and 

E [wr„(0] -> /or a.e. f > 0. (17) 

Finally u satisfies the equation (|TJ) and the boundary conditions @ in the following weak sense: 

/•+oc /• r+co /• r+00 P 

I I d t C,udxdt = / / AMdiv(m(u)DC)dxd£- / / P(u)A(dxdt (18) 
Jo Jn Jo Jn Jo Jn 

for every test function Q G C£° ((0, +00) x Jl) such that D £ • n = on <9f2. 

Remark 1.1. We add a few comments on the previous result: 

• C° ([0,oo); ff 1 (f2)) denotes the space of weakly continuous curves u : [0,oo) — > H 1 ^). 

• A curve u : [0,oo) — > X(Cl) belongs to ACf oc (0, 00; X(Q)) if there exists a function g G 
Lj 2 oc ([0, 00)) such that 



W m (u(s),u(t)) < / g(r)dr for all < s < t < 00. (19) 
• The condition u G AC 2 oc (0, 00; X(fi)) implies that 

< u < M, / u(i, x)dx = m for all t > 0. (20) 



• The condition ([T2"]) implies that the initial datum in ^ is attained in the sense that u(t) 
converges to uo in X(f2) with respect to the distance W m and weakly in if 1 (SI) as t J, 0. 

• Since ue L 2 OC ([0, 00); H 2 (Sl)) and m is LSC, div(m(u)DC) G if oc ([0, +00); L 2 (fl)). 

• Even in the case M = +00, (|15l) . the asymptotic behaviour ([G]), the lower bound (|10l) and 
the Sobolev embedding of H 2 (n) in imply that P(u) G £? oc ([0, +00); L\n)). 

The main examples that fits into the framework of Theorem [T] are the classical Cahn-Hilliard 
equations [I]: For the mobility, one chooses m(r) = r(l — r), so M = 1. Typical choices for the 
free energy G are the double well potential, 

G{r) =dr 2 {l-r) 2 , 

yielding a model for the phase separation for a binary alloy, 

d t u = - div (tt(l -u)D Ait) + 6Au 2 (l - u) 2 ; (21) 

or the function 

G{r) = e(r\ogr + (1 - r)log(l - r)), 6» G R, 
that lead to an equation for the volume fraction of one component in binary gas mixture, 

d t u = -div (it(l - u)D Am) +9Au. (22) 
See [H[5] for a derivation of (f2Tj) and ([22]) . and [14] for a related existence analysis. 
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1.3. The existence result for more general mobilities. The assumption (IM-LSCp about 
mobility can be weakened to condition (M^ 2 I below, at the price that the weak solution to ([T])- 



^ is no longer obtained as a curve of steepest descent in W m , but appears as the weak limit of 
such curves in metrics satisfying (jM-LSCI) . 

Theorem 2. Assume that ft is a smooth bounded convex open subset ofM. d , the mobility function 
m satisfies (|M|) . and the free energy G satisfies In addition, assume that m satisfies 

lims 1/2 m'(s) = 0, and, if M < oo, also limfM - s) 1/2 m'(s) = 0. (M 1/2 ) 

Then, for any initial condition uq £ X(Q) of finite energy E[uo] < oo, there exists a function 
u £ L 2 OC (0, oo; if 2 (ft)) n CJJ, ([0, oo); if 1 (ft)) satisfying the constant mass and maximum estimate 
([2T))l . the energy bound (fT5)l and the weak formulation (fT5|) of equation (JTJ) with the boundary 
conditions ([2]). The initial condition ([3]) is met in the sense that u(t) weakly converges to uq in 
if 1 (ft) as i 4.0. 



The first condition in (M1/2 I is needed to give a meaning to the gradient of m(it) in the 
weak formulation (fT5|) . particularly on the set Z — {(t,x) £ (0, +00) x ft : u(t,x) = 0}. We 
briefly indicate the problem: Since u(t) £ H 2 (fl) for a.e. t > 0, and u(t) satisfies homogeneous 
Neumann boundary conditions, the Lions- Villani-estimate on square roots [22] (see Lemma [A. II 
in the Appendix) states that 

2 <C\\V 2 u(t)\\ L2{n) . (23) 



li 4 (fi) 

Thus u £ L 2 oc (0,oo;_ff 2 (ft)) in combination with dMg/J) implies that m(w) £ i 2 oc (0, 00; ff 1 (ft)) 



In fact, in the proof it turns out that D m(u) = a.e. on the set Z. A similar reasoning applies to 
the zero set of D y/M — u if M < 00. Therefore, it suffices to evaluate the second spatial integral 
in (fTgj) only on the subset {0 < u(t) < M} c ft, thus avoiding to discuss the singularity of m'(s) 
at s — or s — M . Unfortunately, ((23]) cannot be extended to obtain L 2p estimates on roots 
w(£) 1 / p with p > 2, as is easily seen. Without further a priori estimates, there is apparently no 



way to remove condition (M4/2 I- 

The main example in the framework of Theorem [5] is the - (de) stabilized - lubrication or thin 
film equation, where one chooses M = +00 and m(r) = r a with 1/2 < a < 1. The equation is 

d t u = -div(u Q DAw) +kA(/), (24) 

where u : ft x (0,oo) — > [0, +00) describes the height of a thin viscous liquid film on a substrate, 
moving under the influence of surface tension; the lower order perturbation is typically attributed 
to van der Waals forces or similar intermolecular interactions. The destabilized case corresponds 
to k < while for k > 0, the contribution has a stabilizing effect. 

The existence theory of for the unperturbed flow k = is fairly well understood 0. In 
particular, the Hele-Shaw equation obtained for a = 1 has been analyzed thoroughly as a gradient 
flow of the Dirichlet functional in the L 2 -Wasserstein metric, see e.g. [51 [T71 [23] . The perturbed 
flow has essentially been treated in d = 1 dimensions only, see e.g. [21] [28], but some results (e.g. 
on the blow-up behavior of solutions) are available also in multiple dimensions [15j . For non-linear 
mobilities m, the equation's gradient flow structure has apparently not been exploited for rigorous 
analytical treatment before. 

In order to obtain ([M)l from ([T]), one would like to choose 

G(r) = n- r?- a+1 

in the definition of the energy ([5]). This is, however, only possible for certain regimes of /3 and k: 

• If 1 < P < a + 1, then G satisfies (JHJ) for all k £ U. 

• If a < 1 and a + l</3ora = l and /3 > 2 with j5 < 2d/(d - 4) if d > 4, then provided 
that k > 0, i.e., the perturbation must be stabilizing. 

• If P < 1, then there is no way to accommodate the perturbation into our framework. 
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1.4. Key ideas of the existence proof. The discrete approximation scheme in (fTTj) provides 
a family of piecewise constant approximate solutions u T : [0,oo) — > Weak convergence 

towards a limit curve u : [0, oo) — > H 1 ^) along a sequence r„ \. is easily obtained, using the 
machinery developed in [T]. The difficulty lies in identifying the weak limit u as a weak solution 
to ©-©. 

For mobilities satisfying (IM-LSC|) , a semi-discrete version of the weak formulation flT8|) is derived 
by variational methods, i.e., we use suitable perturbations of the minimizers u™ in each step of 
the scheme (|11[) . Our variations of the are obtained by applying an auxiliary gradient flow to 
them. Specifically, in order to arrive at (|18[) , we would like to use variations in the direction of 
the flow generated by the functional 

V[u] := / Vudx 

with a given test function V € C°°(f2) satisfying homogeneous Neumann boundary conditions on 
dQ. To motivate this particular choice, assume for the moment that the test function £ factors as 
C(t,x) = ip(t)V(x). Then the left-hand side of (fT5)l can formally be read as 

/ [ u(t,x)ip'(t)V(x)dxdt = - [ ij(t)^-V[u{t)]dt, 
Jo Jn Jo d * 

i.e., as the temporal derivative of the functional V along the sought gradient flow for E. Further, 

the "flow interchange" Lemma 13.21 taken from [24], says that this expression can equally be 

understood as the temporal derivative of the functional E along the gradient flow of V. Thus, 

variations of the minimizers for (fTTj) along the flow of V are expected to provide a form of (fT5|) . 

Unfortunately, V itself is not a suitable choice for carrying out estimates, since the gradient 

flow generated by V is not regular enough to apply the flow interchange lemma. In particular, the 

functional V is not geodesically A-convex in the metric W m for any A G R (see [7]). As a matter 

of fact, the trajectories of the gradient flow of V with respect to W m are formally given by the 

solutions of the non-linear conservation law 

d t v t = div(m(v t )DV). (25) 

These solutions are expected to develop shocks in finite time. To circumvent this technical problem, 
we consider a modification of V, 

V e [«] := V[u]+eU[u], 
where U is defined in ([7]) , that induces the following viscous regularization in (|25|) 

dtVt = div (m(u t ) D V) + eAvt- 

For LSC mobilities, the viscous regularized flow generated by V e with respect to W m is A £ - 
convex and our strategy goes through. For more general mobilities, even the viscous flow lacks 
convexity. This makes it necessary to perform further approximations: we replace the mobility 
function by LSC-ones, obtain a weak formulation (|18p for the corresponding flows, and then pass 
to the non-LSC limit. 

Even with the discrete version of (|18[> at hand, we still need to facilitate sufficiently strong 
compactness to pass to the time-continuous limit r \, 0. Our key estimate is obtained from the 
dissipation of the functional U in ([7]) along solutions of (JTJ . A direct calculation shows that the 
dissipation term provides a bound in Lj 2 oc (0, oo; H 2 (£l)). The rigorous proof of this if 2 -bound is 
obtained by another application of the strategy above: we interchange flows and calculate the 
variations of E with respect to perturbations of the minimizers in the direction of the gradient 
flow generated by U. This time, the strategy goes through smoothly since the auxiliary functional 
U, which generates the heat flow with homogeneous Neumann boundary condition, is geodesically 
convex with respect to the considered metric W m . 

The paper is organized as follows. Section [5] below provides the essential definitions for the 
measure-theoretic formulation of the problem. In Section [3] we have collected a variety of technical 
results that are applied in subsequent sections to obtain a priori estimates on the semi-discrete 
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approximation u T . Sections |4] and [5] are devoted to the proofs of Theorems [T] and [2 respectively. 
There, we follow the strategy outlined above. 

2. Preliminaries 

2.1. Basic assumptions. Here and in the rest of this paper, we will always assume that 

fl is a convex, smooth and bounded open set of M d , mn := rp-r € (0, M), (fi-conv) 

where M £ (0, +oo] characterizes the domain of the mobility function m as in (|M|). We will always 
assume that m is a mobility function satisfying (|M|) and G is a free energy density satisfying (fGj) . 

2.2. Notation: admissible and regular densities. As in ©, for a given mobility m and a 
mass m satisfying (|f2-conv|) we introduce the sets of admissible and regular densities on f2 

X(Q) := jpG-L^fi) 0<p<M, Jpdx = mj, (26) 

X r (ft) := {p e C°°(f7) < inf p < supp < M, Jpdx = Xtt^. (27) 

Since we will keep fixed the mobility m (and its domain of definition (0,M)) and the total mass 
m, we will omit to indicate the explicit dependence of the above spaces from these two parameters. 

We often identify an element u £ X(Cl) with the nonnegative measure u — u££ d in M. d supported 
in f2 and we will consider weak convergence of sequences in X(Cl) in the sense of distributions of 
&'{M. d ): 

u n ^u in %>'(K. d ) lim / u n <pdx= [ uipdx for all <p G C™(R d ). (28) 



n— >oo 



Remark 2.1. Since u n are nonnegative with fixed total mass, we could also equivalently consider 
test functions (p £ C°(M d ) in (f2"5| ; when M < oo, X(£l) is a convex and bounded subset of L°°(f2) 
and (|28[) also coincides with the weak-* convergence in the latter space. 

For every extended-valued real functional £F : X(Cl) — > (— oo, +oo] we denote by Dom(5") its 
proper domain Dom(3 r ) :— {u £ X(£l) : ^(u) < oo}. J is called proper if Dom(3 r ) is not empty. 

We will consider curves in X r (£l): they are maps 7 : [0, 1] — > X(il) which we will also identify 
with functions 7 : [0, 1] x Q -> [0, M] such that j(t, ■) e X r (fi) for every t £ [0, 1]. We say that 

7 : [0, 1] -> X r (f2) is regular if f ^ 7 (t, x) £ ^([0, 1]) for all x £ Q. (29) 

In a similar way, a functional 

J : X r (f2) -> K is regular if i ^ J[ 7 (i)] e C 1 ([0, 1]) for every regular curve 7, (30) 

and a map S : [0, 00) x X r (il) — > X r (fl) is regular if the curves 

S(-,it), S(t, 7 (-)) are regular for every u £ X r (£l), t>0, and for every regular curve 7. (31) 

2.3. Survey: weighted transport distances. We shall now review the weighted transport 
distances W m introduced in [T2] (see also [7] and [53]) without going into details about their 
formal definition. When M — 00 they could in fact be pseudo-metrics, i.e. they satisfy all the 
axioms of the usual notion of distance except for the fact that the value +00 may be attained; 
nevertheless, even in the case M = 00 the next proposition shows that the restriction of W m to 
the sublevels of the convex functional (recall (0) 

U[u] := [ U{u) dx where U"{s) = — ^, U{m n ) = U'{m n ) = 0, (32) 

is a finite distance. Notice moreover that, besides m, W m also depends on the domain f2: we will 
denote it by W mi o when we want to stress this dependence. In particular, for every 8 > we will 
also sometimes consider the 5- neighborhood Qui of f2 

n [s] -.= n + m d = {x e R d \ dist(x, n) < 5}. (33) 
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Proposition 2.2. In the setting of §2.1[ the pseudo-metric W m on the space X(Q) has the 

following properties: 

(a) For every u ni u G X(f2), 

lim W m (n n ,!i)=0 <=> u n ^-u in&{M. d ) (34) 



according to (|28[) f&wi see afeo remark [Kl 
(b ) For every c > the sublevels of U 



|u G X(f2) | U[u] < c| are compact metric spaces w.r.t. W m . (35) 

(cj For every decreasing sequence of convex sets Q" converging to £1, if two sequences Uq, u" G 
converge to uq and u% in the sense of distributions respectively, i.e. 

lim / v?<pdx= [ uupdx for all (p G C?°(M. d ), i = 0, 1, (36) 
™^°° Jo" 7n 

i/ien 

W m ,n(«o,«i) < liminfW m , "K\<). (37) 

n— 7-+oo 

(g^ 7/7 : [0,1] — > X r (0) is a regular curve according to (|29[) and y> : [0,1] — > is i/ie 

corresponding curve of weak solutions to 

— div (111(7) D tp) — d s j in Q, n-Dip — on <9f2, (38) 

i/ien i/ie W m -distance between uq = 7(0) and u\ = 7(1) is bounded as follows: 



W m (u Q , Ul ) 2 < f /m(7)|D^| 2 dxds. 
Jo Jn 



(39) 



lim 

n— ► oo 



fej Let ?ii G X(fi), i = 0, 1, oe given with W m (uo,ui) < 00. For every decreasing sequence of 
smooth convex sets f2„ converging to fi as ?i -> 00, such that fl n D &[8 n ] f or a vanishing 
sequence 5 n , there exists a sequence of regular curves ("approximate geodesies") 7„ : [0,1] — > 
X r (Q n ), such that: 

• 7n(0) and 7«(1) converge to uq and U\, respectively, in L 1 (R d ) as n — > 00 and for every 
proper and lower semicontinuous convex integrand F : [0, M] — > [0, 00] 

F(j n (i,x))dx= / F(ui)dx i = 0,l. (40) 
Jn 

• if (p n : [0,1] — > TJ^fin) are the corresponding curves of weak solutions to (|38[) on Q, n , 
then 

W m (u ,wi) 2 = lim / / m( 7 „)|D^„| 2 da;ds. (41) 

No properties of the metric W m other than those listed above will be used in the sequel. Notice 
that (|59"|) &; (j4"l"j) establish the connection between the metric W m and the formal definition of the 
metric tensor given in ([§])■ 

2.4. The entropy and energy functionals. In this section, we derive some relevant properties 
of the entropy and the energy densities U, G introduced in ([7]) and ©. For definiteness, we make 
the following specific choice for the function U in ((?]): 

f s s — r m 

"<*> "/„=«*' "» : = m " = M- (42) 

Lemma 2.3. The entropy functional U is lower semi-continuous with respect to the weak conver- 
gence (|28p , and satisfies the following bounds 

< U[u] < C(l + ||u||| 2(n) ), for all u G Dom(U). (43) 

The constant C above only depends on fl, m, and m(so) = m(mn). 
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Proof. Lower semi-continuity is a consequence of the convexity of U, which, in turn, follows from 
U"(s) = l/m(s) > for every s € (0, M). The lower bound in (14"5|) follows from non- negativity 
of U, indeed convexity of U and yield that mo is a minimum for U and U(mn) = 0. For 
showing the upper bound in (|4*5)) . first note that 



*(so) 



s if s < sq 



^m(so) if s > so, M = +oo, 

by concavity of m. Thus for G := and G M := ^(s") 

^(aQogs - log(s ) - 1) + s ) < G if < s < s , 

f/(s) < \ s^((Af - s)(log(M - s) - log(M - s ) - 1) + M - s ) < C M if a > *o, M < ™, 
2^)( s ~ s o) 2 < 2^) s2 if s>s , Af = oo. 

Now (|4^]) follows by the boundedness of f2. □ 

Concerning the function G, we decompose its second derivative G" into the difference of its 
positive and negative part 

L :=G" = L+ - L-, L~ := — min(G",0), L+ := max(G", 0). (45) 

Fixing s € (0, M) (e.g. so = mn as before) and assuming without loss of generality that G(s ) = 
G'(so) = (recall that the integral of elements in X(Cl) is fixed to be m) we have the decomposition 

G = G conv + G conc , G conv (s)= f L+(r)(s-r)dr, G conc (s) = - f L~(r)(« - r) dr, (46) 

and the corresponding one 

P = Pmc, + Pdcc,, P mci (s)= f L+ {r)m{r) dr, P d ecr(s) = - / L~(r)ni(r) dr. (47) 

(gBJ) and the upper bound L~ < G(l + 1/m) = G(l + U") of ^ yield 

G C onv(s) > G conv (s ) = = G conc (s ) > G conc (s) > -C(l + s 2 ) for all s G (0, M), (48) 

proving the lower bound in ([TU1) . It follows immediately from the lower bounds in JG| that Pdecr 
is Lipschitz continuous and G con c is continuous in [0,M) (and also in M if M < oo) since it is 
concave and bounded from below. 

In order to check the continuity of G conv in (the same argument applies to M when M < oo), 
let us first observe that P[ ncT — L + m is integrable around since Pi ncr = P — Pdccr is locally 
bounded around by ([G]). Recalling (PHI) we easily get for < s < so 

r s ° P! (r) sn f So 

G co „v( S ) = / -^Y(r -s)dr< / P/ ncr (r) dr < ~P incr (0). 

Js m W m(s ) 7 S 

Since G con v is convex we conclude that it has a right limit at 0. 

With (flu| and the above remarks at our disposal, we can obtain simple lower bounds on the 
the energy functional E defined in ©. 

Lemma 2.4 (Basic properties of E). The functional E is bounded from below in the space X(Q) 
and lower semi- continuous with respect to the distributional convergence (|28|) in the space X(ii). 
Finally the following estimate holds 

1 



[ G conv {u) dx < E[u] + E for all u e X(Cl), (49) 
Jo. 

where G conv has been defined by (|46[) - (|45p and the constant Eo only depend on Q, the mass m, the 
dimension d and the function G. 
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Proof. To begin with, we prove the estimate (|49[). We recall the following Gagliardo-Nirenberg 
( 16 . 25 ) interpolation inequality 

NU 3 (n)<C , i||D W ||i 2(n) |H|^ n) +C' 2 ||ix|| il(n) , ueH\n) (50) 

where 9 = d/(d + 2) and the constants C\, G2 only depend on f2 and d. In our specific case of 
u £ X(O) n ^(O) we have, for every e > 0, 

||«IU-(n) < Ci|| D uWi^xa 1 - 6 + G 2 m < e\\ T>u\\ L 2 {n) + G 3 (e)m (51) 

where G 3 (e) := (2C 1 /e) 1 /( 1 ~ e '> + G 2 . In particular 

||«||| a( n) <2 £ 2 ||D w || 2 2(0) +2G 3 ( £ ) 2 m 2 . (52) 

Using the decomposition (14*6")) . the lower bound (|38"|). and (|52l . for the constant G in flU^) we have 



E[«] > h Du|| 2 2 - C(|fi| + || U || 2 2 ) + / G conv (u)dx 
z Jn 

> (i -2e 2 c)||D w || 2 2 -G(|fi| + 2G 3 ( £ ) 2 m 2 )+^G conv (u)d ; r 



(53) 



Choosing e 2 = 1/(8C) in ([53]) and using again ([52]) with e 2 = 1/2 we obtain (gH]) with the 
constant E := C(\Q\ + 2G 3 (l/(2v / 2G)) 2 m 2 ) + l/4G 3 (l/\/2) 2 m 2 . Boundedness of E from below 
is an immediate consequence of (|49p . recalling that G conv is non- negative. 

In order to prove lower semi-continuity, assume that a sequence Uk £ X(fl) converges to a 
limit u £ X(fl) according to (|28|) . It is not restrictive to assume that u k £ H 1 ^) and that 
sup fe ^ M E[«t] < +00. By estimate (|49|) the sequence u k is bounded in _ff 1 (J7). Hence, up to 
subsequences, Uk converges weakly in H 1 ^), converges strongly in L 2 (tt), and converges pointwise 
_Sf d -a.e. to u. 

By (|4"5)l and Fatou's Lemma we have that 

liminf/ G conc (u k ) + G(l + u 2 k )dx > / G conc (u) + G(l + u 2 )dx. (54) 



The L 2 (n) strong convergence of Uk and concavity of G con c yield 



limsup / G conc (u k ) + C(l + u%)dx < / G conc (u) + G(l + u 2 )dx. (55) 
fc->+oo Jn Jn 

From JM]) and (JSS]) it follows 

lim / G conc (uk )dx = / G conc (u)dx. (56) 
k^+coJn Jn 

Second, by Fatou's Lemma, it follows that 



liminf / G conv (u k )dx > / G conv (u)dx. (57) 

fe^ + oo J Q J n 



Finally, since Duk converges weakly in L 2 (il) to Dm, 



k —>■-(- 00 



n Jn 



liminf / |Du fc | 2 dir> / |Du| 2 dx. (58) 



The lower semi-continuity of E follows from (|5"6")) , (|57p and . □ 

We will denote by E m ; n the minimum value (depending on f2, m, G) of E on X(fl). Notice that 
estimate (|4"3"|) in combination with (j4"9"]l yields 

< U[u] < G(E[w] + E ), for all it € Dom(E), (59) 

with some constant C only depending on m(so), fi, m and G. 
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3. A PRIORI ESTIMATES 

3.1. Semi-discrete approximation. We begin by invoking a result from [1 that guarantees the 
well-posedness of the minimizing movement scheme (jlll) . i.e. the existence of the semi-discrete 
curves u T and their compactness for vanishing step size r 4- 0. 

Proposition 3.1. In the setting of H2.ll for every uo G X(fl) n Dom(E) and r > there exists a 
sequence {u"}„>o satisfying and the following energy estimate: 

1 N 

E[u^] + — J2^ m K,u^ 1 ) 2 <-E[u ] forallNeN. (60) 

n=l 

Moreover, for every sequence r„ 4- i/iere exists a subsequence, still denoted by r n , and a continuous 
limit curve u : [0, +oo) — > X(Q) such that U Tn (t) converges weakly to Ut in for every t > 0. 

The curve u is globally 1/2-Holder continuous 

W m («t,u.) < (2(E[ Uo ]-E min )) 1/2 | s -i| 1 / 2 /or aZi s,t6 [0, +oo). (61) 

The curve tt-*u t satisfies ([15]). ([T5|). (fT71) and (|2Dj) . 

Proposition 13.11 is obtained by combining the results from Chapters 2 and 3 (see in particular 
Sections 2.1, 2.2 and Corollary 3.3.4) of Q]. The properties of E proven in Lemma [2^41 are sufficient 
to apply this general theory and the uniform upper bound on E[u T ] given by (|60p improves the 
pointwise convergence of U Tn with respect to W m to the weak convergence in H l (£l). It should 
be remarked that we do not claim uniqueness of solutions, even on this discrete level, except in 
the case when E is a convex functional. 

3.2. Flow interchange lemma. For the derivation of r-independent a priori estimates on the 
interpolations u T , we employ the device of the flow interchange lemma, which has been proven in 
[24] . Before reviewing the lemma and its proof, we recall the definition of A-flow in the metric 
space X(SY) given in [8]. 

Definition 1. Let SF : X(Q) — > (— oo,+oo] be a proper lower semi- continuous functional and 
X G K. A continuous semi-group S* : Dom(If) — > Dom(3 r ), t > 0, is a A-flow for 2f if it satisfies 
the Evolution Variational Inequality (EVI) 



— lim sup 

2 /ij,o 



W m (S h (u),v) 2 - W m (u,v) 2 



h 

for all measures u,v G 001x1(5") with W m (u, v) < +oo. 



^W m (u,v) 2 + J[u] < J[v], (62) 



Recall that a continuous semigroup S on a set D C X(Q) is a family of maps S* : D — > D, 
t > 0, satisfying 

S t+S (u) = S*(S s (u)), limW m (S t (u),u) = for all u € D. (63) 

Notice that the continuity of S is already coded in (j()2"j) : it is sufficient to choose v := u in (|r?2"j) . 
(f63|) and the triangle inequality yields 

W m (S t (u),w) < +oo for aU u G D, t> 0; (64) 

in particular the "lim sup" in (|6"2"|) is well defined. 

Lemma 3.2 (Flow interchange Lemma |24|). Assume that is a X-flow for the proper, lower 
semi- continuous functional 5F in X(Q) and let u™ be a n-th step approximation constructed by the 
minimizing movement scheme (|11[) . If u™ G Dom(9 r ) then 

3> T ] - J[u T ] < rhmmf I I - -W m (w r ,u T ) . (65) 
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Proof. (C.f. [24]) By definition of it™ as a minimizer in (fTT|) . 

^W m «,<- 1 ) 2 + EK] < i-W m (S^«),<- 1 ) 2 + E[S^«)] 
holds for every h > 0. This implies 
-rliminf [hT 1 (E[S£«)] - E[<])] < -Urn sup [ft- 1 W m (S t ? «), u^ 1 ) 2 - W m «, u?" 1 ) 2 ] . 

To conclude (|S5|) from here, apply (j()2"j) with the choices u = u™ and v = □ 



3.3. Eulerian calculus. In order to apply the flow interchange Lemma 13.21 with a particular 
auxiliary functional £F, we need to exhibit the associated semigroup S (usually given implicitly 
as the solution to a nonlinear evolution equation) and to verify that it is indeed a A-flow, i.e., it 
satisfies the EVI (j62|) with a finite constant A. A very general strategy to attack this problem 
is the Eulerian calculus for transportation metrics, that has been developed by the third author 
in [5], based on earlier work by Otto and Westdickenberg [57] • Similar to the flow interchange 
estimate, the basic idea is to simplify estimates by exchanging two time-like derivatives. We also 
need that S can be suitably approximated by semigroups on smooth densities: here is the relevant 
definition. 

Definition 2. Let us fix a nonnegative vanishing sequence 5 nt let Q n D ^[5 ra ]> a, decreasing sequence 
of smooth convex sets converging to Q, let J : X(Q) — > (-co, +oo] be proper and l.s.c. functionals, 
and let S be a semi-group on Dom(3 r ) C X(Q,). 

We say that {3>i, S ra } ne jsj is a family of mollifications for 2r, S if 

(a) 3>j : X r (n„) — » M is a regular functional according to (|30[) . 

(b) Each S„ is a regular semi-group on X r (Q n ) according to (|31l) . 

(c) For every po, pi € Dom(5") with W m (/3o,pi) < oo the regular densities 7n(i) given as in (e) 
of Proposition \2.2\ satisfy 

lim J„[ 7 „(i)] = 3[pi] i = 0,l. (66) 

(d) For every sequence p n £ Dom(J n ) converging to p € Dom(3 r ) in L 1 (R d ) as n — >■ oo with 
3n[p«] — >• ^[p] we have S^p„ S*p in Sl'iW 1 ) and liminfV^oo ? n [S^p„] > U^S'/o] for all 
t > 0. 

Let 7, 3 n and S, S„ as in the definition above. For a given n £ N consider a regular curve 
7„ : [0, 1] — ► X r (Sl n ) and for every h > 0, introduce 7* : [0, 1] — > X r (yi n ) — the h-perturbation of 
In — by 

Notice that the (h,s) i-> 7^ l (s, x) is of class C 1 in [0,oo) x [0,1] thanks to the regularity of S n . 
Also, introduce the action of the perturbed curves 

A h n (s)= [ |D^( s )| 2 m( 7 £( S ))dz, (67) 

where the V?n( s ) ^ ^ 1 (^n) form a s-differentiable family of solutions to the associated Neumann 
problems 

-div(m(^(«,-))D^(»,.))=^75(*,-) mfi„, n-D^( s ,-)=0 on 9Q„. (68) 

These Neumann problems are solvable because 7^ is a regular curve of densities in X r (£!„); in 
particular, the mass is constant, and thus 9 s 7„(s, •) has vanishing average on fl n . The following 
result is essentially an adaptation of Theorem 2.2 in [5] to the situation at hand. 

Lemma 3.3. Under the hypotheses and with the definitions above, assume that inf„ 3 n > F > —00 
and h i-» 3>j[S„p] are non-increasing for every p € X r (H, n ), and there exists A < (independent 
of n and of the considered curves j n ) such that the inequality 

ld h A h n (s) + s\A h n ( s ) < -d s J n [^(s)} (69) 
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holds for all s G [0, 1] and all h > 0. Then S is a X-flow for 3". 



Proof. The core idea is to prove and integrated form of (|62[) by estimating the perturbed action 
(p7|) starting from a family of approximating geodesies - provided by (jgj in Proposition 12.21 - 
connecting two given admissible measures u, v G Dom(5") c X(£i). Without loss of generality, we 
assume that 5F n are non-negative and A < 0; the case A = follows by obvious modifications. 

Given u, v G Dom(iF) C X(Q) at finite distance, and a family of approximating geodesic 7„ on 
X r (f2„) between v and u in the sense of (jej) in Proposition [221 define w° = 7„(0) and u\ — 7 n (l). 
Multiply (|69j) by e 2A,l;5 and integrate with respect to s G [0, 1]; this gives 



l_d_ 

2&h 



f e 2Xhs A h n (s)ds < - f e 2Xh *d s ? n [-yZ(s)]ds 
Jo Jo 



2\he 2Xhs ? n [ri(s)]ds 



= J n [u° n ] - e 2M J„[S^] + / 

Jo 

< J n [<]-e 2 ^J„[S£<], 
since A < while 5F„ is non-negative. Next, integrate with respect to h G [0, if], which yields 

- Jf e 2 ^A«( S )d.s + — — ?„[SM<- y A°„( S )d S + J ff^[<], 

where we also used the fact that h i— > 3>j[S„u*] is non- increasing. Further, a reparametrization of 
s H> 7„(s) similar to that used in [8j Lemma 5.2] yields 

-2XH _ i 

1 -W m «,S^) 2 < / e 2 ^A^( s )d S . 
Jo 



-2Af? 
In summary, we have 



,-2Aif _ i I _ 2XH i /-i 

2Air W m «,SK) 2 + -^-J n [SM<- y A°„( s )d S + J ffJ„[<]. 

By our choice of 7„, (e) of Proposition 12.21 and (ISfJl) yield on one hand that 

lim / A°„(s)ds = W m (u,v) 2 , lim ? n [<] = %]. 

On the other hand, we know by Proposition 12.21 and the properties listed in Definition [2] that it* 
converges to the density u in L 1 (M d ) and 

3{S H u] < liminf ? n [S?ul\, W m {S H u,v) 2 < liminf W m (Sf ul , u° n ) 2 . 
Altogether, this yields the inequality 

-2XH _ 1 1 _ 2Aff i 

— ^-W m (S ff u, v) 2 + _ 6 2X ?lS H u] < g W m («, «) 2 + H3H 
from which the EVI property (fr?2")) is deduced after division by if > in the limit if I 0. □ 



4. Proof of Theorem Q] 
Throughout this section we use the notation introduced in ^ 1 . 2L 
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4.1. iJ 2 -regularity and strong convergence. The goal of the following is to prove: 

Proposition 4.1. In the setting of each solution u™ of the minimizing movement scheme 
([lip satisfies 

< £ H 2 {Cl) for all n € N, r > 0, (70) 
and the following bound holds for the piecewise constant interpolant u T , 
T 

\\u T (t)\\ 2 H2(n) dt < CT(E Q + ~E[u Q }) <+oo forallT>0, (71) 

with a constant C independent of t > 0. Moreover, every sequence \. contains a subsequence 
(still denoted by Tk) such that 

u Tk -> u strongly in L 2 (0, T; H 1 ^)) for all T > 0, (72) 

u Tk -> u weaA% m L 2 (0, T; i? 2 (0)) /or a/Z T > 0. (73) 

The proof of Proposition 14.11 rests on the fact that the densities u T are r-uniformly bounded 
in L 2 (0, T; H 2 (fl)) for arbitrary T > 0. As motivation for the arguments below, we provide the 
relevant formal calculations in the case G = 0: assuming that u is a smooth solution to ((T|) 
satisfying (|2j). differentiation of the entropy functional U[u(i)] introduced in Q yields 

— U[u(t)]= / U'(u(t))d t u(t)dx = - I U'(u(t))div(m(u(t))BAu(t))dx 
d* Jo Jn 

= / LT"(u(i))m(ti(i))Du(t)-DAu(i)do;= / B u(t) ■ B Au(t)dx 
Jn Jn 

= - [ (Au(t)fdx 
Jn 

because of the identity U"(s) — l/m(s). Consequently, U[u(t)] is decreasing with respect to t, 
and (still formally), 

U[u(T)]+ / / (Au) 2 da;dt < U[« ]. 
Jo Jn 

Taking also into account the contribution of G and the convexity of CI, one ends up with an estimate 
of the form (f7Tj) . The goal for the rest of this section is the rigorous proof of this estimate. 

We wish to apply the flow interchange Lemma |3~21 with = U. To this end, we need to identify 
the associated semi-group S, with S t v given by the smooth solution v t to the Neumann problem 

dtVt = Avt in Ct, n • Bv t — on dft, vo — v. (74) 

Lemma 4.2. The semi-group S induced by solutions Vt of the problem (|74[) on X r (fl) extends to 
a 0-flow S for U. 

This fact is a special case of a more general result proven in Theorem 6.1]. We provide the 
relevant calculations for the specific situation of Lemma B~2l as we shall refer to it later. 

Proof. We wish to apply Lemma l3.3l In order to define a family of mollifiers {J n , S„} ne N for U, S, 
we consider a sequence of domains fl n := f2[i n i for some vanishing sequence S n > and we define 
3>j(u) := U„[u] := J„ U (u(x))dx and the heat-semigroup S„ on fl n with homogeneous Neumann 
boundary conditions. The contraction properties of the heat flow imply that convergence of the 
initial conditions in L 1 imply the same convergence of the solution at any time h > and it is 
easy to verify all properties required in Definition [5] 

Let regular curves j n : [0, 1] — > X r (£l n ) be given. By classical parabolic theory, the /i-perturbed 
curves 7^ : [0, 1] — > X r (fl n ) are well-defined for any h > 0, and for every s £ [0, 1], the function 
(h, x) h-> Jn(s, x) is a classical solution to 

d h ^ = sA^ inO n , D 7 £-n = on dfl n , 7 °(s) = 7n (s). (75) 
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We need to verify the principal estimate (|69p of Eulerian calculus, which reads in the situation at 
hand (we will omit to indicate n in the following) as follows: 

~ / a /l [m( 7 ' 1 )|D^| 2 ] dx < - [ d s [U( 7 h )] Ax. (76) 
^ jo Jn 

Using the definition of ip h in (|68[) and its boundary conditions, the right-hand side evaluates after 
integration by parts to 

- f U'{-/ h )d s ^ h dx = - [ U"( 1 h )D 1 h - (m{'j h )Dip h )dx = - f Dj h -D<p h dx = f y h Aip h dx 
Jn Jn Jn Jn 

since U" (s)m(s) = 1 by definition of U. For the /i-derivative of the action, we find 

\ f d h [m( 7 h )\B V h \ 2 ] dx = \ [ d h m( 7 h )\T><p h \ 2 dx + [ m( 7 h )d h D <p h ■ D ip h dx 

A Jn z Jn Jn ^\ 

= -~ / 9 h m( 7 ' l )|D^| 2 dx+ / (d h m( 1 h )Di P h + m( 1 h )d h D<p h ) -Dp" da:. 
* Jn Jn 

To simplify the second integral above, first observe that for every smooth function 9 € C°°(f2) it 
follows from that 

f m{j h )D<p h -D9dx= [ d sl h 9dx. (78) 
Jn Jn 

Taking the /i-derivative yields 

f (d h m(~f h )D6 + m{i h )d h D(p h )-'Dddx = [ d h d s j h 6dx. (79) 
Jn Jn 

If 9 satisfies homogeneous Neumann boundary conditions, n • D 9 = on <9S1, we obtain from (|75l) 

d h -yHdx = s I -) h Mdx, (80) 
! Jn 

and the s-derivative amounts, in view of ((75)1 . to 

/ d s d h y h <&dx= j ^Atidx + s [ m(j h ) D <p h ■ D Atf dx, (81) 

JO JO JO 

and thus allows to express the mixed derivative dhd s ^ h in (|79p . Using as test function 6* = ip h in 
((79")) and (|5T)l . the integrals in (1771) become 



i f S A [m(7*)|D/| a ] te = -\l d h m( 7 h )\D<p h \ 2 dx 

+ / /A^'dx + s / m^D^-DAt/dx. 
JO jo 

We evaluate the first integral on the right-hand side, 

-\ f d h m(<y h )\D<p h \ !, dx = -± f m , (7 h )«A7 fc |D^| a da; 
1 Jn 1 Jn 

= S -[ D(m'( 7 ,l )| D(f h \ 2 ) ■ D~f h dx 
1 Jn 

f D (m'( 7 ' l ) N ) ■D"/ h \Dip h \ 2 dx + - f m'^ h ) D / • D (| D <p h \ 2 ) dx 
Jn ' 2 J n 

f m"( 7 /l )|D 7 'f|D^Tdx- J / m( 7 /l )A(|D^| 2 )d2: 

JO 2 Jf2 

I / m( 7 h )D(|D^f)-nd<^ d 
z jao 



(82) 



in * Jn 

'so 

Using the last identity in (|82|). taking into account the Bochner formula 
— -A(| D ty9 /l | 2 ) + DAi/) 11 • Dtp' 1 = — 1| D 2 /|| 2 < 0, 
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and that D (| D (p h \ 2 ) ■ n < on <9f2 since f2 is convex, see (|124[) . we find that 

~ / d h [m( 7 h )\-Dip h \ 2 ]dx<- f d s [U(-f h )]dx+ S - f m'V^Dyf |D^fd.T. (83) 
1 Jq Jn z Jn 

By concavity of m, this proves ([75)1 . □ 

The following Lemma provides the last missing piece for proving (|71l) by means of the flow 
interchange Lemma F3.2[ namely the dissipation of the energy E along the heat flow (ITU) . 

Lemma 4.3. Let v : [0, oo) -> H 1 ^) be a solution to (f73|). J/ 

hminf -(E[v s ] - E[v }) > -oo, (84) 

siO S 

then vq € H 2 (fl) and 

-hminf -(EK1-EW) > \ I (Av f dx - C(E + E[v }), (85) 
where the constant C depends only on m(so) , |0| and G. 

Proof. By classical parabolic theory, the solution v to ([74]) is smooth, and for every < so < 
si < oo it satisfies < inf x£ f2 v s (x) < sup xen v s (x) < M for (x,s) € O x (so,Si). Thus E[w s ] is 
continuously differentiable with respect to s in [so, si] with 

-r-E[u s ] = [ Dv s -DAv s dx+ f G'(v g )Av s dx=- f (Av s ) 2 dx - / G"(v s )\ Bv s \ 2 dx, (86) 
ds Jn Jn Jn Jn 

where the last equality follows after integration by parts, using that the boundary condition 

n • Dw s = is satisfied for any s > 0. Taking into account ([G]) the second integral can be 

estimated as follows, 

- / G> s )|D« s | 2 dx<- / G'l(v s )\Bv s \ 2 dx<c( [ \Dv s \ 2 dx+ [ l -^^dx 
Jn Jn Wn Jn m(u a ) 

Recall dH]) and the identity 

v s 1 1 

In case that M = +oo, we obtain 



\Bv s \ 2 J ^ s f \Bv, 12 



dx < 7 ; / J — -^-dxH ^— / |DuJ 2 da; 



m(v s ) m(s ) Jo «s m(s ) 

< c(^|DV^| 2 dx + ^|D Us | 2 dx). 
Moreover, by Holder's inequality and estimate (11 23[) from the Appendix, 

f | d ^ri 2 dx < f i d v^ri 4 a*) 1/2 < ( (d+ |j )|n| / (a« s ) 2 dx) 1/2 < e f (a Vs ) 2 dx + k e . 

with X e = (d + 8)|f2|/(64e). And analogously, if M < oo, 

D V ° 1 2 dx < « I M d , + / | D(M — t; s )| 2 ^ 



m(u a ) m(s ) J n w s ' m(s ) M - u s 

<C( / |D>I| 2 dx+ / |Dv/M-i; s | 2 dx), 



where we use that 



/ \D^/M -v s \ 2 dx<e j {Av s ) 2 dx + K e 
Jn Jn 

Choosing e above sufficiently small, and observing that 

/ |D« s | 2 da; < / |Dw | 2 da;, 
Jn Jn 
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-EN < ~J^Av s ) 2 dx + C{l + |KH^ ( n)) 5 



it thus can be achieved that 
_d 
d? 

for a suitable constant G. 

Recall that the curve s i— > v „ is continuous in if 1 (SI) and that G can be decomposed as in 
The continuity of G conc and the lower bound yield that 



lim / G conc (w g ) dx = / G conc (v ) dx; 
<4° Jn Jn 



in Jn 
on the other hand, since G conv is convex, we have 



Gconc(w s ) d.x < / G conv (v )dx for all s > 
n Jo 



so that Fatou's Lemma and the continuity of G yields 

lim / G conv (v s )dx = / G CO nv(vo)dx. 

Consequently, the function s t— > E[u s ] is continuous at s = and we have that 

i(E[w,]-E[«o]) <-~ / (A^ (s) ) 2 d.T + G(l + || Wo || 2 Hl(f ,)) ) 

with < 9(s) < s. By it follows that the family {Av e ^ s - ) } se ^ >So ) for so > is weakly compact 
in L 2 (S1). Since w s converges to vq strongly in as s I 0, we have that vo £ H 2 (Vl) and 

-hminfi(EH-EN) > ihminf f (Av e(s) ) 2 dx - C(l + \\v f H . {Q) ) 

> ^(A^cLr-^l + IHI^n)). 
Another application of the estimate finally provides (|85D . □ 



Proof of Proposition \4-l\ By Lemma 12.31 we can apply the flow interchange Lemma 13.21 with £F = 
U. By Lemma T4.3I applied to vq — it™ we have that u" lies in H 2 (£l), and by and (|65p . for 
any nePj, it follows that 

J / (A<) 2 dx < UK" 1 ] - U[<] + G(E + E[«?])t. (87) 

Here the constant G is the same as in (|55|). and does not depend on r, on n or on the solution u T . 

Let T > and r e (0, 1) be given, and define N G N such that (N — l)r < T < Nr. In view 
of (|59p and E[u™] < E[uq], summing (|87p from n — 1 to n — N, we find that the interpolating 
function u T satisfies 

W^r\\h {0 . T]LHn)) <rJ2 / (A<) 2 dx<2GT(E + EK]), 
n=l ^° 



which is obviously independent of t € (0, 1). Combining this with (|^|) and applying again (|122p , 
we conclude that u T remains uniformly bounded in L 2 (0, T; iJ 2 (SI)) as r | 0, for any T > 0: 



/ IM*)llfl*(n) d * ^ ^(Eo + E M) < +CX3. 

Jo 



By (|88[) we have that, up to subsequences, w Tre converge weakly to u in I/ 2 (0, T; if 2 (Sl)) for every 
T > 0. Since we have already seen in Proposition 13.11 that u T pointwise converge weakly in 
1 (SI) and thus strongly in L 2 {VL) by Rellich's Theorem, the dominated convergence theorem 
shows that u T converges strongly in L 2 (0, T; L 2 (S1)). From here, strong convergence u Tn —> u in 
L 2 (0, T; iJ x (S7)) follows by standard interpolation between the uniform bound (|88p. □ 
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Corollary 4.4. In the setting of Q2.1\ we have for all T > 

P{u Tn ) -> P{u) strongly in L^O, T; L 1 ^)), (89) 

and, if m is also Lipschitz ( as for (|M-LSC[) ) 

m(u T „) -» m(u) strongly in L 2 (0, T; H 1 ^)). (90) 

Proof. is trivial when M < oo. When M = +oo, by Sobolev imbedding Theorem and the 
uniform estimates (|7ip and (|S(J)) we know that 

J J [[u T f + \u T \)dxdt<C T (91) 

uniformly with respect to r. Since u T (up to subsequence) converges strongly to u in i 1 ((0, T)xfl), 
we deduce the same property for P(u) thanks to (|G|) . 

(|M| is a standard consequence of the fact that m is Lipschitz and C 1 . □ 



4.2. Weak formulation. The remaining section is devoted to prove the following Proposition 
stating that the time-continuous limit u obtained before is a weak solution in the sense of (|18[) . 

Proposition 4.5. Under the assumptions of Theorem^ let V be a spatial test function satisfying 

VeC°°(n), BV-n = on 8Q } (92) 
and a temporal test function tj) £ C^°(0, +oo) be given. Then 

r+oo r+oo 

i/)'(t)V[u(t)]dt= ij(t)N[u(t),V]dt, (93) 
Jo Jo 

where the nonlinear functional N is given by 

N[u,V] :=- Audiv(m(u)DV)dx+ / P(u) AVdx. 
Jn Jn 

In the spirit of the ideas developed in |19] , we would like to use the flow interchange Lemma 
13.21 with J := V the potential energy functional V : X(fi) — > K defined by 

V[«] := / V(x) u(x) dx, 
Jn 

with a test function V satisfying ((92|) . As already mentioned in the introduction, the functionals V 
are — unfortunately — never A-convex (for any AGl) along geodesies of the space (X(fl), W m ), 
unless the mobility m is a linear function |7, Section 2.3]. To cure this problem, we shall construct 
a A e -flow for the regularized functional 

V e [«] := V[u] +eV[u], (94) 

with e > instead, which amounts to solutions of the classical viscous approximation of ([231) , 

d s v s - div (m(u s ) D V) - eAv s = in (0, +oo) x fl with Dv s ■ n = on (0, +oo) x dn. (95) 

Proposition 4.6. Under the assumptions (|M]h (jM-LSCp on m, suppose that V satisfies ([9^| . 
Define the semigroup S e by taking S|vo — v s , the unique solution to (|95[) with initial condition vq. 
Then S £ extends to a X E -flow S e for V e with respect to W m , with some X £ > —K/e where K > 
only depends on V and m. 

Proof. As in the proof of Lemma 14.21 we need to verify (|59"]) for the flow S e and the functional 
£F = The calculations are similar to the proof there, but more terms need to be controlled. 

Below, we shall implicitly use various properties of the solution semi-group S £ for f|95|) . A 
summary of these relevant properties are given in Lemma IA.2I in the Appendix. In particular, 
note that S e is well-defined and L 1 -continuous on the admissible densities X(Q), and that it 
leaves the regular densities A r (fi) invariant. 

In order to define a family of mollifications {V e n , S e ,„} for V £ , S e , we assume without restriction 
that G f2, we take a monotone sequence r) n i 1 and we define fl n := rj n Cl = {rj n x : x £ O}, 
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choosing S n I so that fl[s n ] C f2 n . We define V^(x) = V(x/rj n ) and V Ein [u] = / n V^(a;)da; + 
eU[u]. Then V n satisfies (|9"2f in f2 n We define for every n the solution semi- group S £-n of the 
problem l|95|> for V n on the domain fi n . It is not difficult to check that all the conditions of 
Definition [5] are satisfied. 

We turn to prove (|69p . writing for simplicity il in place of Q n everywhere. The s-derivative of 
V e amounts to 

- / d s [eU(-y h ) + V-y h ] dx = e [ j h A(p h dx + [ m( 7 /l ) D V ■ D <p h dx. 
Jn " Jo Jn 

Moreover, the weak formulation ([50)1 is modified as follows, 

/ d hl h $dx = se I -f h Atfdx-s j m( 7 /l ) DV ■ Dtfdx, (96) 
Jo Jn Jn 



and, consequently, (I81|) is replaced by 

/ d s d hl h $dx = e [ y h A#dx + se [ m( 7 ' l )D/ • D Atfdx 
Jn Jn Jn 

- / m(-f k )DV -Dtfdx-s / m(-f k )D(p k ■ D (m' (y h ) ~D V -Dtf). 
Jo Jo 

Performing the same manipulations as in the proof of Lemma 14.21 one obtains 

1 / 9 fe m( 7 ' l )|D^| 2 da ; = ^ / m"( 7 h )|D 7 fc | 2 |D/| 2 dx-se / m( 7 ' l )A(| D 7 'f) da; 

2 jo 2 jq Jn 

-s / m( 7 ' i )m'( 7 ' i )Dl/-DV l -D( y 5 /i da:-s f m{ 1 h )m"( 1 h ) D 7 ' 1 • D V\ D ^| 2 dx. 
Jo Jo 

Summing up everything provides 

~ [ 5 h [m(V , )|D^| 2 ] dx = -c> s ( / V^dx + e [ U(j h )dx) 
2 Jo V Jo Jo / 

+ y^m"(/)|D 7 'f |D^| 2 dx (97) 

+ | / m( 7 A )D(|D/| a ).ndJr«'- 1 (98) 
^ Jao 

-s [ m( 7 ' 1 )m"( 7 ' i )(D 7 ' l -D^ ,l D^ ,l -Dl/-D 7 ' l -Dt/|D^ i | 2 )dx (99) 
Jo 

-s I m{j h )m'(j h )Dp h D 2 VF>v h dx. (100) 
Jo 

We need to show that the sum of the terms from (|!JT)) to (llOOp are less than — sX £ A for a sufficiently 
small (negative) constant A e . The integral in (I100[) is readily controlled by a multiple of A, recalling 
that m has the Lipschitz property (|M-LSC[) and observing that 

|m( 7 ' l )m'( 7 ' l )D^ l D 2 ^D^| < sup |m'( S )||| V||c=(0) m( 7 ' l )|D^f • 

s 

In order to absorb the integral in (|M)l into the (non-positive) integral in (|97p and a multiple of A, 
we apply Young's inequality to the integrand and estimate 

|m( 7 ' l )m"( 7 ' l )(D 7 ' 1 ■ D Lp h D Lp h ■ DV — D -f h ■ DV\D ip h \ 2 )\ < 2m( 7 /l )|m"( 7 ' l )|| D-y fc || D <p h \ 2 \ D V\ 
<|(-m'V))|D7 h | 2 |D/| 2 + -(-mV)m( 7 ,t ) 2 )|D^| 2 |Dy| 2 



<|(-m''( 7 h ))|D 7 ^| 2 |D^| 2 + ^|| 2 7l(0) sup(-m''( S )m( S )) m(/)|D^| 2 . 



Thus, defining, for every n 



2 

A e ,„ := -sup|m'(s)|||y„||c=(n„) - -||^n|lc?i(n„) sup(-m"(s)m(s)), 

s £ s 
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and recalling that (1981) is non-positive for convexity of fi we obtain 

If 9, l [m( 7 ' l )|D^| 2 ]d.T + S A £ ,„ / [m^^lD^'fjdx 
z Jq„ Ju„ 



< -a 



(101) 



Defining A e := inf„ A e . n > -co (thanks to the uniform boundedness of all the derivatives of V n ) 
the principal estimate follows from (|101l) . □ 

The flow interchange estimate (|55|) is applicable. To obtain a sensible a priori estimate, we still 
need to express the dissipation term in (|65[) . 

Lemma 4.7. Let v s be as in Proposition \4-6\ and assume that vq € Dom(E) (1 H 2 (Q). Then 

- liminf 1 (E[«J - E[«q]) > -N[v 0l V] + e(i J^o) 2 dx - C(E + E[«o])) . (102) 

Proof. For 5 > sufficiently small, define approximations of G by 

' G (s+M^a) HM<oo, 
^G{5 + s) ifA/ = +oo. 



G s (s) 



This regularizes the possible singularities of G'(s) for s J. and s f Af. Denote by Ea the energy 
functional associated to Gs instead of G. 

For the following calculations, we need some properties of solutions to the problem ([95]) . which 
are summarized in Lemma lA. 21 in the Appendix. By (I125P we have that s M- Ea[u s ] is absolutely 
continuous and we can calculate for almost every s > its derivative 

-E s [v s ] =e(- J Av s Av s dx + J G' s (v s )Av s dx] 



ds 



< -e 



[ Av s div(m(u s ) D V) dx + [ G' s (v s ) div(m(« s ) D V) dx 
Jn Jn 

\J {Av s f dx - C(E + E tf [«o])) + N s [v s , V}. 



,2 

The last estimate is obtained by treating the term multiplied by e exactly as in the proof of 
Lemma 14.31 and integrating by parts in the last two integrals (which is allowed for the smooth 
approximation Gs and does not produce boundary terms since V satisfies homogeneous Neumann 
conditions (|9"2"|) V Moreover, following the proof of Lemma |4~5I is is easy to check that the constant 
C in the last integral can be chosen uniformly with respect to 5. Then we have 

EJuJ - E s [v ] ^ 1 f ( e 



-\j ( ~ l X (Awt)2 dx + Ns[vu dt + eC(i + E5[w ° ]) ' 



By (JTUJ) it is easy to check that J Q G s (v)dx -> J n G(v)dx and J n P s (v)AV dx -> j n P(v)AVdx 
as 8 I 0. Passing to the limit as S j. we obtain 

EH - E H ^ 1 [* ( e 

s 

By the right continuity property (|126[) we can pass to the limit by s \. obtaining (| 102[) . □ 

The flow interchange estimate ([55)1 provides the following. 

Lemma 4.8. Let V be a given test function satisfying (1921) and ip £ C£°(0, +00) be a given 
temporal test function satisfying ip>0. Then, 

^'(t)V e [u T (t)]dt< / t/j T (t)N[u T (t),V]dt + C-, (103) 



where the simple function ip T : (0, +00) — » [0, +00) is defined by ip T (t) — ip((n — l)r) for (n— l)r < 
t < nr for all n G N. The constant C in (|103[) is independent of t and e and depends only on the 
test functions V and ip, and on the initial energy E[uq]. 
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Proof. Since u T is a simple function with respect to t > 0, which is constant on intervals ((n — 
l)r, tit], and ^ is smooth with compact support, it follows that, for some sufficiently large N £ N, 

- / v'(0v £ [« r (t)] dt = - V / ^(t)v e [<] dt 

JO n=1 J(n-l)r 
N 

= X;W(n-l)r)-V(nr))V B [u?] (104) 

n=l 

AT 

= ^^((n~l)r)(V £ K]-V e [ U r 1 ])- 

n=l 

By Proposition 14.61 we can apply the flow interchange Lemma 13.21 with 5F = V e . By (|7D|) we can 
apply Lemma 14.71 with vq = u". Combining the flow interchange estimate (1651) and inequality 
(|TU2)| we find 

VeKJ-VeK- 1 ] <rN[< ) y] + £ r((7(Eo+E[uo])-iy'(A<) 2 dx)+^W m (<,<- 1 ) 2 . (105) 
Combining (| 104|) with f| 105[> and recalling that i/i>0we obtain 

,+oo N N-l 



/ i>'{t)V e [u r (t)} dt<rY, i>{{n - l)r)N[u?, V] + eC(E + E[u ])r £ ^(nr) 

n=l n=0 
+ 00 

p^(t)r^ (- 



_l supr-'- 1 

n— 1 

/■+°° _ /•+<» 

</ ^ T (t)N[u T (t),V r ]di + C(E +E[ Mo ])e / ip T {t)dt 
Jo Jo 

Kt 

+ sup^(t)(E[u ]+Eo), 

where the energy inequality (joT)]) has been used to obtain the last line. The claim (| 1031) follows. □ 

In order to finish the proof of (|93l) , we pass to the time-continuous limit r 1 and the limit as 
e \. simultaneously, in such a way that the remainder term in (|103D goes to zero. 



Proof of Proposition \4-5\ For definiteness, let (r n ) n ^ be a vanishing sequence for which u Tn — > u 
strongly in L 2 (0, T, iJ 1 (il)) according with Proposition 14.11 Without loss of generality, for (|72|) 
we may further assume that u Tn — > it and Dtl Tii -)D« almost everywhere on (0, +oo) x f2. It is 
sufficient to choose the vanishing sequence e n := -^/tvT in order to have that Cj 1 - 1 in (1103[) . 

We start by proving convergence of the left-hand side in (|103|) . By the bounds from ([59]) and 
the the monotonicity of the energy (jBT)]) . one finds that 

\V s [u T {t)]-V[u T {t)]\ < eC(E + E[uo]) (106) 

for every r > and t > 0. Choosing T > such that supp(-0) C [0, T], using (|106l) . we have 

ip'(t)V en [u Tn (t)]dt- [ <i/>'(t)V[u(t)]dt 



< sup W{t)\\ I \V £n [u Tn (t)}-V[u Tn {t)}\dt+ I \V[u Tn (t)]-V[u(t)]\dt 

te(o,T) 







dt + 


Jo 



< sup |V/(t)| (Te„<7(Eo + E[ti ]) + sup|V(a:)| / 
t£(o,T) y xeo Jo Jn 



u| dx dt 



which shows that 

r+co r+co 

lim / if}'(t)V en [u Tn (t)]dt= / ^'(t)V[u(t)] dt. (107) 
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From (|103p and (|107[) one concludes that 

p+co r+co 

ip'(t)V[u(t)}dt < liminf / tp T (t)N[u T (i), V] dt. (108) 
Jo ™^°° Jo 

Next, we claim that the minimum limit in (|108[) is actually a limit, and that 

lim / $ T Au T div (m(u r ) D V) dx dt = \ ipAu div (m(u) D V) dx dt, (109) 

"->+<*> J sir " Jo r 

lim / 4) T P l {u T )AVdxdt= / ^(^AVdadi, (110) 



for i = 1,2. In fact, (I109p and (IllOp follow almost immediately from Corollary 14.41 Combining 
(1901) with the weak convergence (|73[) and the uniform convergence of tp Tn to tp in ^t, one obtains 
(|109l) . And recalling that tp Tn uniformly converges to ip in Ot, we obtain (II 10[) from (|8"9"|) . 
Inserting (ITM| and ([TTUl) into (fTUSj) we obtain that 

x/}'(t)V[u(t)]dt< ip(t)N[u(t), V] dt (111) 



for all V G C°°(f2) satisfying (|92l) . and all non-negative ip G C£°(0,+oo). Exchanging V with 
— V in (|lllj) yields the respective equality (|93p. Trivially, extends from non- negative test 
functions ip to all ip G C%°(Q, +oo), thus hnishing the proof. □ 

Since any space-temporal test function in (fT8|) ( G C°°((0, +oo) x f2), with D ( ■ n = on dil, 
can be approximated in C°°((0, +oo) x f2) by sums of functions of the type ((t,x) = ip(t)V(x) 
with V G C°°(IT) satisfying and V> G C c °°(0, +00), Theorem □follows. 

5. Proof of Theorem [2] 

5.1. Approximation. Theorem[5]is now proven by approximation of the more general mobility 
function m satisfying (M4/2 I by mobilities that have the Lipschitz property (|M-LSCp . To this 
end, define for all 8 > sufficiently small: 

• if M < +00: 



mis) :=m(^i S + 4)-^ 



where Sg < s 2 s are the two solutions of m(s) = 6. 
• if M = +00: 

mg{s) := m(s + sj) - S. 

where s$ > is the unique solution of m(s) = S. 
Introduce accordingly Pg by 



P s (s) = / ms(r)G"(r)dr. (112) 



Lemma 5.1. For all 6 > sufficiently small, the mj are smooth functions that have the Lipschitz 
property (|M-LSCp and satisfy the pointwise bounds < < m. In particular, we have 

W m (u,«) < W m5 (u,w) for allu,v e X(n). (113) 

-For <5 J, 0, t/ie mj converge monotonically and globally uniformly to m. Moreover, if G satisfies 
(iGl) wii/i respect to m, f/ien it afeo satisfies (|G|) wii/i respect to each mj. Finally, the P$ are 
continuous functions, and there is a constant K such that 

-K(l + s 2 ) < P 5 (s) < P(s) + K(l + s) (114) 

for all s G (0, M) and all S > sufficiently small, and P$ converges to P as 5 \. 0, uniformly on 
[0, M] if M < 00, or uniformly on each [0,s] if M = +00. 
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Proof. Smoothness, non- negativity and the Lipschitz property of mj are evident from its definition, 
and the concavity and smoothness of m. In the case M = +00, also the upper bound m«5 < m is 
a trivial consequence of concavity, as is the uniform convergence for 6^,0: 

< m(s) - m s (s) <S- m'(s + s s )s s < S. (115) 

In the case M < +00, the upper bound can be proven as follows: assume that m attains its 
maximal value at a £ (0, M); then m$ attains its maximum at as = {a — s\)M / '(s 2 — Sg). Without 
loss of generality, assume as < a. For all s £ [0, as], we have 



<(s) = ^j/* m' ( "l^s + s\)< m'00 



>s 

and thus also ttig(s) < m(s). A similar argument provides 1115(5) < m(s) for all s £ [a,M]. For 
s £ [era , a] , the inequality follows since is non-increasing and m is non-decreasing on that 
interval. The argument for uniform convergence of 111,5 to m is established essentially with the 
same argument as in (|115[) . making again a case distinction whether s £ [0,(75], s £ [ss,a], or 
s £ [a,M\. 

Condition fG} on G is less stringent for the approximations since mj < m in case M = +00, 
and 1115/(1 + m<s) < m/(l + m) in case M < 00. Concerning the continuity of Ps, we remark 
that, in view of m d - < m, the integrability of m(s)G" (s) near s = (and near s = M if M < 00) 
implies the respective integrability of m.g(s)G"(s). Moreover, 

sup \P(s) - P s (s)\ < [ (xn(s)-m s (s))\G"(s)\ds 

0<s<s Jo 

in combination with the pointwise convergence of 1115 to m implies uniform convergence of Ps 
to P on all intervals [0,s]; notice that the dominated convergence theorem is applicable since 
< (m — m<s) |G" | < m|G"|, and the latter is integrable by assumption. Finally, if M — +00, 
then (| 1 14[) is another consequence of ([G]). Indeed, on one hand, 



Ps(s)= m 5 (r)G"(r)dr > -C / (1 + m(r))dr > -K{1 + s 2 ), 
Jo Jo 

and on the other hand, using also (|115[) . 

Pg(s) = P(s) - f (m(r) - m s (r))G"{r)dr < P(s)+S f (G"(r))_dr 
^0 Jo 

< P(s) +CS [ ( 1 + -^r-r)dr < P(s) + K(l + s). 
Jo v m(r)/ 



For M < 00, (I114p simply amounts to <5-uniform boundedness of Ps, which is clear from the 
uniform convergence to P. □ 



5.2. Weak and strong convergence. Lemma 15.11 implies that Theorem [T] is applicable to 
the approximate mobilities for each i5 sufficiently small: there exist respective solutions 
us : [0, +00) -> X(Q) of problem i.e., 

(•+00 p+oo /• P+OO p 

I / d t (u s dxdt= / / Au s B(m s (us)BC)dxdt+ / / Ps{u s )A(;dx dt (116) 
Jo Jn Jo Jn Jo Jn 

for all test functions £ £ C£°((0, +00) x il) such that D £ • n = on <9f2. We wish to pass to the 
limit as S I in ([TT6t . 

Lemma 5.2. There exists an absolutely continuous curve u : [0,oo) — > X(fl) and a vanishing 
sequence Sk such that the us k converge to u weakly in ff 1 (Sl) and strongly in L 2 (Q) pointwise in 
time, as well as strongly in L 2 (0,T; L 2 (fl)) for every T > and pointwise a.e. on (0, +00) x il. 
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Proof. The curves us satisfy estimate (|6lj) in the respective metric W ms with the global Holder 
constant determined by E[uo]. In view of (|113l) . the family (us)s>o is equi-continuous with respect 
to W m . Moreover, the sublevel sets of the energy E are compact. The claim on convergence now 
follows by the variant of the Arzela-Ascoli theorem given in [I] Proposition 3.3.1]. An application 
of the dominated convergence theorem with respect to time provides strong convergence of us 
to u also in sense of L 2 (0, T; L 2 (£l)), and thus (without loss of generality) also pointwise a.e. 
convergence. □ 

In the following we write 5 J. to indicate "along a suitable vanishing sequence 5k for k — > oo" . 
Lemma 5.3. For every T > 0, 

limsup / \\us(t)\\ 2 H 2 (n) dt < +oo. (117) 
sio Jo 

Consequently, us converges to u weakly in L 2 (0, T; H 2 (fl)) and strongly L 2 (0, T; H (Cl)) as 6 J, 0. 
Moreover, Ps(us) converges to P(u) in L l (Q, T; L (Q,)). 

Proof. Define the ^-approximations of the entropy functional Us as in (|4"i21) with m^ instead of 
m. Then estimate ((59)) holds with a constant C independent of S for every Us, at least for 
all S > sufficiently small. Indeed, observe that m^so) > ^m(so) if 6 is small enough, and 
hence inequality (|44l) follows. In the same way, inequalities (1551) an d (EFJ can be obtained with 
^-independent constants C. In combination, (|117l) follows. 

The stated weak convergence of us is now a consequence of Alaoglu's theorem and the uniqueness 
of the weak limit. The strong convergence is obtained by interpolation of the strong convergence 
in L 2 (0, T; L 2 (S1)) with the bound (|117l) . To prove convergence of Ps(us), we argue as in Corollary 

mi □ 

5.3. Convergence of the mobility gradient. The next goal is to establish convergence of 
m s (u s ) in L 2 (0,T;H l (Cl)). 

Lemma 5.4. ms(us) converges to m(u) strongly in L 2 (Q,t) for every T > as 6 4- 0. 

Proof. The uniform convergence of the mobility functions mj to m, and the pointwise a.e. conver- 
gence of us to u suffice to conclude pointwise a.e. convergence of ms(us) to m(u) on fix- Moreover, 
if M < oo, then ms(us) is (5-independently bounded, and by dominated convergence it follows 
that ms(us) converges strongly to m(u) in L 2 (yix). In the case M = +oo, the ^-uniform bound 
ms(s) < m(s) < C(l + s) and the strong convergence of us in L 2 {Qt) imply equi-integrablity of 
Ini^M^)! 2 in Q,t- We invoke Vitali's theorem to conclude the proof. □ 

For the proof of convergence of the gradients Dmjfiij) in L 2 (VLt), we distinguish the cases 
M < oo and M = oo. 



Lemma 5.5. Assume M < oo. Define g : [0, M\ — > R by gs(s) := \J s(M — s)m' s (s) for all S > 
sufficiently small. For S \. 0, the gs converge uniformly to the continuous function go : [0, M] — > M 
given by go(s) = v/ s(M — s)m'(s) for all s G (0, M), and go(0) — go(M) = 0. 



Proof. Let e > be given. Since \J s{M — s)m'(s) —> for s I and for s f M, respectively, by 



hypothesis (M 1 / 2 ), and since m'(s) is an non-increasing function on (0,M), there are a\ < a 



2 



such that \gs(s)\ < \J s{M - s)\m'{s)\ < e for all s e [0, a\] U [a 2 , M] , and all 5 > sufficiently 
small. For S 4- 0, the points s$ and s 2 in the definition of mj converge to and to M, respectively. 
By smoothness of m, one thus has local uniform convergence of m^ to m', and consequently 
\gs{s) — go(s)\ < e for all s e [o-\ , a 2 ], for sufficiently small S > 0. The uniform convergence of the 
gs also proves continuity of go- D 



Lemma 5.6. Assume M < oo. Then Drnj^) converges to Dm(u) strongly in L 2 (Q,t) for every 
T > as S I 0. 
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Proof. Observe that 

\Btn 5 (u 5 {t))\ 2 = ( gs ( Us )) 2 -L(\ D V^WI 2 + I D y/M - u s {t)\ 2 ) (118) 

for every t > at which ug(t) £ H 2 (Q). Pointwise convergence of us to u almost everywhere on 
fix and uniform convergence of gs to go imply pointwise convergence of the compositions gg(ug) 
to go(u) almost everywhere. In combination with the <5-uniform boundedness of gs it follows in 
particular that gs(u$) — > go(u) in L 4 (f2*r). 

Now let Q C fix be a measurable subset. From (II 18[) we obtain 



/ | Dmg(ug) | 2 dx dt = / |g A -(u 5 )| 2 | D ^/ua| 2 dx di + -!- / |g 5 (-« 5 )| 2 | D a/M - us\ 2 dxdt 
Jq m Jq m Jq 

- ^(/ T |i D ^wn^ d *) ' {j Q m{ - us)i dxdt ) 1 (119) 

1/2 * ' ' 

|| D «*Wllt* dtj g s (u s fdxdt 



1 

+ M 

Choose Q = Zt := {(t, a;) G f2<r : u(t,x) — or u(t,x) = M}. Since gs(us) — > go(u) = in 
L 4 (Zt), the right-hand side of (|119[) vanishes as 8 1 0, and Dmj(«j) — >• in L 2 (Zt). 

Further, since the Q in inequality (II 19|) can be chosen as any arbitrary subset of fix \ Zt, the 
equi-integrability of \gs(us)\ 4 on Qt is inherited to |Dm«5(u5)| 2 on {It \ Zt- On Ot \ ^t, the 
composition m'(u) is everywhere finite, and m' s (us) — > m'(u) pointwise a.e. In combination with 
the pointwise a.e. convergence of Dua, it follows that 

D ms(us) — m' s (us) D —> m'(u) Du = D m(u) 

strongly in L 2 ({It \ Zt), invoking Vitali's theorem once again. 

The proof of strong convergence mg(ug) — > m(u) in L 2 (0, T; H 1 (fi)) is now concluded by 
observing the following. Since mg(ug) — > m(it) in L 2 ({It), the strong L 2 (f2x)-limit of Dms(us) 
coincides with the (uniquely determined) distributional derivative Dm(w). In particular, Dm(u) 
vanishes a.e. on the set Zt- □ 

Lemma 5.7. Assume M = +oo. Then Dmjfuj) converges to Dm(w) strongly in L 2 ({It) for 
every T > as S 10. 

Proof. Similar as in the proof of Lemma [5.61 define gs : |0,oo) — > K by gs{s) := y/sm' s (s) for 
all 5 > small enough. Analogously to Lemma 15.51 one obtains convergence — uniformly on 
every interval [0; s] — of gs to the continuous function go : [0, oo) — > R with go(s) '■= v / * m '( s ) f° r 
s € (0, +oo), and go(0) = 0. Like in (I118p . we observe that 



\-Dms(us(t))\ 2 = 4(gs(u s )) 2 \D y/u^l 2 (120) 
for every t > at which us(t) € H 2 (fl). As before, we conclude that gs(us) converges to go(u) 



almost everywhere on (0, +oo) x Q. Moreover, by construction of the m^ and (M^l, one has 
< gs{s) < C(l + s 1 / 2 ) for all s > 0. Since \us\ 2 is equi-integrable in Qt for arbitrary T > 0, 
the compositions 1 3,5 (71,5) | 4 are also equi-integrable in {It- By Vitali's Theorem, it follows that 
gs(us) -t go{u ) in L 4 (VL T ). 

From (|120|) . we conclude that 

/ \Dm.s(us)\ 2 dxdt < A i \gg(ug) | 2 1 D ^/us | 2 dx dt 
Jq Jq 

<*(jf II D V^WIIi* d *) ' (/ gs(u s ) 4 dxd?j ' (121) 



holds for any measurable set Q C {It- From this point on, the proof is identical to the one for 
Lemma with the only change that Zt ■= {(t, x) € {It ■ u(t, x) — 0}. □ 
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5.4. End of the proof. We are now able to pass to the limit in the weak formulation (11 16|) . By 
Lemma [??2l Aus converges to Au weakly in L 2 (Ht) as 6 I 0, and by Lemma l5~6"l or l5.7l respectively, 
Dms(us) converges to Dm(w) strongly in that space. Hence, the term inside the first integral on 
the right-hand side of (|116l) converges weakly to the desired limit in Convergence of the 

second integral follows from the strong convergence of Ps(us) to P(u) in L}(SIt) stated in Lemma 
15.31 This finishes the proof of Theorem [2j 

Appendix A. Appendix 

We recall a Sobolev like inequality, that will be useful in order to estimate the rate of dissipation 
of the energy E along the heat flow. 

Lemma A.l. Assume that fl is a smooth convex open set and that u 6 H 2 (fl) satisfies homoge- 
neous Neumann boundary conditions, n ■ D u = on dfl. Then 

[ ||D 2 u|| 2 dx< f (Au) 2 dx<d[ ||D 2 u|| 2 da;. (122) 
Jn Jn Jn 

If, in addition, u is non-negative, then y/u £ W 1,4 (fl) and 

16 / | Dv^| 4 dx < (d + 8) / (Au) 2 dx. (123) 
Jn Jn 



Proof. By density, it obviously suffices to prove (|122p for u £ C°°(Sl) satisfying n- Du = on dil. 
Integrating by parts, it follows 

D 2 w|| 2 dx= / n-D 2 7i-Diidcr— / D Au -Dudx 
Jdn 



/ n-D 2 u-Dudo--/ n-DuAudcr+ / (Au) 2 dx. 
Jan Jan Jn 



'an Jon 
Now observe that the second boundary integral vanishes because of the no-flux boundary condi- 
tions, whereas the integrand in the first is pointwise non-negative, 

n-D 2 w-Du<0. (124) 

(For a proof of this classical fact, see e.g. |18l Lemma 5.2].) This shows the first inequality in 
(11221) . The second inequality follows by observing the pointwise relation 

d d d 

(Au) 2 = duud jjU < - {d u u) 2 + (d n u) 2 = dY, (duu) 2 < d\\ D 2 u\\ 2 . 

i,j=l i,j = l i=l 

The estimate (|123[) follows by combination of [TSJ Lemma 3.1] with (I122[) above. □ 



The next lemma summarizes selected properties of solutions to the viscous conservation law (I95[) . 
These are used in the calculation of the right derivative of the energy E along the corresponding 
flow. 

Lemma A. 2. Assume that il is a smooth bounded domain, m satisfies ([Mjl. (|M-LSCp and V 

satisfies (|92j) . If vq € X r (Q), then there exists a unique smooth classical solution v(s,x) of 
problem (|95l) and v(s, ■) € X r (tt). If vq and vq are two initial conditions in X r (Q), then the 
L 1 -contractivity estimate 

IK - "slUi(n) < ||«o - uo||ii(fi) 

holds at any time s > for the associated solutions v, v. If vq £ H 2 (Q), then v s £ H 2 (Q) for 
every s > 0, 

Av, div(m(u) D V) £ L°°(0, +oo; L 2 (n)) } (125) 

and the maps 

s i y Av s , s i — y div(m(w s ) D V) are right continuous with values in L 2 (Q). (126) 
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Proof. This lemma is deduced with standard methods for parabolic equations, see e.g. |20) . and 
we leave most of its proof to the interested reader. Here, we shall only comment on the well- 
definiteness of the flow on X r (fl). From the classical theory, it follows that v G C°°(IR+ x fl), i.e., 
the generalized solution is smooth and classical for positive times, and is continuous initially. One 
then easily verifies, using (jM-LSCp . that v satisfies the comparison principle. Next, observe that 
(I95p admits a family of stationary solutions (w c ) c£ k of the form 

v c ( x ) = F _1 (c-e _x V(x)) with F(s) = 

It is easily seen from the properties of m and by smoothness of V, that v c € X r (fl) for every 
c G R, and that v c \. or v c f M uniformly on SI, for c j -oo or c | oo, respectively. Hence, 
any solution v can be sandwiched between two stationary solutions. Since vq g X r (Q), one has 
v(s) G X r (fl) for all s > 0. 

□ 
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